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Abstract 

We reconsider the homogeneous Faddeev-Merkuriev integral equations for three- 
body Coulombic systems with attractive Coulomb interactions and point out that 
the resonant solutions are contaminated with spurious resonances. The spurious so- 
lutions are related to the splitting of the attractive Coulomb potential into short- 
and long-range parts, which is inherent in the approach, but arbitrary to some ex- 
tent. By varying the parameters of the splitting the spurious solutions can easily be 
ruled out. We solve the integral equations by using the Coulomb-Sturmian separa- 
ble expansion approach. This solution method provides an exact description of the 
threshold phenomena. We have found several new S-wave resonances in the e~e + e~ 
system in the vicinity of thresholds. 
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1 Introduction 



Certainly, to calculate resonances in three body atomic systems the methods 
based on complex coordinates are the most popular. The complex rotation 
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of the coordinates turns the resonant behavior of the wave function into a 
bound-state-like asymptotic behavior. Then, standard bound-state methods 
become applicable also for calculating resonances. The complex rotation of 
the coordinates does not change the discrete spectrum, the branch cut, cor- 
responding to scattering states, however, is rotated down onto the complex 
energy plane, and as a consequence, resonant states from the unphysical sheet 
become accessible. 

In practice bound state variational approaches are used, that results in a 
discretization of the rotated continuum. By changing the rotation angle the 
points corresponding to the continuum move, while those corresponding to dis- 
crete states, like bound and resonant states, stay. This way one can determine 
resonance parameters. 

However, in practical calculations the points of the discretized continuum scat- 
ter around the rotated-down straight line. So, especially around thresholds it 
is not easy to decide whether a point is a resonance point or it belongs to the 
rotated continuum. 

Recently, we have developed a method for calculating resonances in three-body 
Coulombic systems by solving homogeneous Faddeev-Merkuriev integral equa- 
tions [1] using the Coulomb-Sturmian separable expansion approach [2]. As a 
test case, we calculated the resonances of the negative positronium ion. We 
found all the 12 resonances presented in Ref. [3] and observed good agreements 
in all cases. 

With our method we succeeded in locating 10 more resonances in the same en- 
ergy region, all of them are very close to the thresholds. This is certainly due to 
the fact that in Ref. [2] the threshold behaviors are exactly taken into account. 
Unexpectedly we also observed that in the case of attractive Coulomb inter- 
actions the Faddeev-Merkuriev integral equations produce numerous resonant 
solutions of dubious origin. We tend to regard them as spurious resonances 
which come about from the somewhat arbitrary splitting of the potential in 
the three-body configuration space into short-range and long-range terms, that 
is an inherent attribute of the theory. 

Since the possible appearance of spurious states in the formalism is new and 
surprising, for the sake of the better understanding of its mechanism, in Sec- 
tion 2 we outline briefly the Faddeev-Merkuriev integral equation formalism. 
In our particular case of the e~e + e~ system two particles are identical and 
the set of Faddeev-Merkuriev integral equations can be reduced to an one- 
component equation. In Section 3 we discuss the spurious solutions. In Sec- 
tion 4 the Coulomb-Sturmian separable expansion method is applied to the 
one-component integral equation. In Section 5 we present our new calculations 
for the L = resonances of the e~e + e~ system and compare them with the 



2 



results of the complex scaling calculations of Ref . [3] . 



2 Faddeev-Merkuriev integral equations 

The Hamiltonian of a three-body Coulombic system is given by 

H = H° + v° + v° + v°, (1) 

where H° is the three-body kinetic energy operator and v% denotes the Coulomb 
potential in the subsystem a, with a = 1,2,3. We use throughout the usual 
configuration-space Jacobi coordinates x a and y a . Thus v%, the potential be- 
tween particles f3 and 7, depends on x a . 

The Hamiltonian (1) is defined in the three-body Hilbert space. The three- 
body kinetic energy, when the center-of-mass motion is separated, is given 
by ' 

H° = hl+hl = hl + hl = hX + hl, (2) 



where h° is the two-body kinetic energy. The two-body potential operators 
are formally embedded in the three-body Hilbert space v = v c (x)l y , where 
l y is a unit operator in the two-body Hilbert space associated with the y 
coordinate. 

Merkuriev proposed [1] to split the Coulomb interaction in the three-body 
configuration space into short- and long-range terms 

vZ = vP + vW (3) 



where the short- and long-range parts are defined via a splitting function: 



v ( a s) =v^C{x a ,y a ) 



v° t 1 - C(x a ,y a ) 



(4) 
(5) 



The splitting function ( is defined in such a way that 



KmJ(x,y) = { 



1, if \x\ <x°(l + \y\/yy/», 
0, otherwise, 



(6) 
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Fig. 1. Typical shape of for attractive Coulomb potentials, 
where x°,y° > and v > 2. Usually the functional form 

C(x, y) = 2/{l + exp [(x/x°Y/(l + y/y )} } , (7) 

is used. So, the separation into short- and long-range parts is made along a 
parabola-like curve over the {x, y} plane. Typical shapes for v ^ and are 
shown in Figures 1 and 2, respectively. 

In atomic three-particle systems the sign of the charge of two particles are 
always identical. Let us denote in our e~e + e~ system the two electrons by 
1 and 2, and the positron by 3. In this case v% , the interaction between the 
two electrons, is a repulsive Coulomb potential which does not support two- 
body bound states. Therefore the entire v% can be considered as long-range 
potential. With splitting (3) the Hamiltonian can formally be written in a 
form which looks like an usual three-body Hamiltonian with two short-range 
potentials 

H = H^+v{ s) + v ( 2 s \ (8) 

where the long-range Hamiltonian is defined as 

H®=H° + v?+v$ ) +v°. (9) 
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Fig. 2. Typical shape of for attractive Coulomb potentials. 

Then, the Faddeev method is applicable and, in this particular case, results 
in a splitting of the wave function into two components 

l*) = liM + hfe>, (io) 

with components defined by 

|^> = G%)^|*>, (11) 

where a — 1, 2 and G^(z) = (z — if z is a complex number. 

In the case of bound and resonant states the wave-function components satisfy 
the homogeneous two-component Faddeev-Merkuriev integral equations 

\^) = G^{z)v?\^) (12) 
\^) = Gl\z)vi s) \^) (13) 

at real and complex energies, respectively. Here G^} is the resolvent of the 
channel long-ranged Hamiltonian H$ = + G®(z) — (z — if®) -1 . 
Merkuriev has proved that Eqs. (12-13) possess compact kernels, and this 
property remains valid also for complex energies z = E — iT/2, T > 0. 



Further simplification can be achieved if we take into account that the two 
electrons, particles f and 2, are identical and indistinguishable. Then, the 
Faddeev components j^i) and \ip2), m their own natural Jacobi coordinates, 
have the same functional forms 

(xml^i) = (x2y 2 \ifa) = (xy\ip). (14) 

On the other hand 

W=P%), (15) 

where V is the operator for the permutation of indexes 1 and 2 and p — ±1 
are eigenvalues of V. Therefore we can determine from the first equation 
only 

|^) =G?vi 8) pV\rl>). (16) 

It should be emphasized, that so far we did not make any approximation, 
and although this integral equation has only one component, yet it gives full 
account both of asymptotic and symmetry properties of the system. 



3 Spurious resonance solutions 

Let us suppose for a moment that the y° parameter in C( x iiVi) i s infinite. 
Then ( would not depend on y±, the separation of the potential into short 
and long range parts would go along a straight line, and v[ l \xi,yi) would be 
like a valley in the y\ direction. The potential vf\xi, yi), which is Coulomb- 
like along xi, would support infinitely many bound states and, as free motion 
along the coordinate y\ would be possible, the Hamiltonian H° + would 
have infinitely many two-body channels. 

If, however, y° is finite, the straight line along the y± direction becomes a 
parabola-like curve and the valley, as y± goes to infinity, gets broader and 
broader and shallower and shallower, and finally disappears (see Fig 2). As the 
valley gets closed the continuum of H° + associated with the y\ coordinate 
becomes discretized. So, if y° is finite, H° + have infinitely many bound 
states. Similar analysis is valid also for H° + v%\ and consequently, also 
has infinitely many bound states. 

This, however, due to (11), can lead to spurious solutions. If, at some energy, 
has bound state has pole. Then, \ip a ) in Eq. (11), irrespective of 



6 



can be finite even if \^/) is infinitesimal. These solutions are called spurious 
solutions: although the Faddeev components are not identically zero, but their 
sum vanishes. 

Let us examine now the spectral properties of the Hamiltonian 

H® = H® + v{ s) =H° + + v$> + v%. (17) 

The three-body potential is attractive and constructed so that (^2,2/2) 
vanishes as y 2 tends to infinity. Therefore, there are no two-body channels 
associated with fragmentation 2, and of course neither with fragmentation 3, 
the Hamiltonian Hi has only 1-type two-body asymptotic channels. 

What happens to the bound states associated with v^l They are embedded 
in the the continuum of H° + v f , and become resonant states. So, by solving 
the Faddeev-Merkuriev integral equations, at some complex energies, we can 
encounter spurious solutions. These spurious solutions are not related to the 
Hamiltonian H, but rather only to some auxiliary potentials coming from the 
splitting procedure. Consequently, the spurious resonances should be sensitive 
to the parameters of (, while the physical resonances not. This way we can 
easily distinguish between physical and spurious resonance solutions. 



4 Coulomb-Sturmian potential separable expansion approach 



We solve Eq. (16) by using the Coulomb-Sturmian separable expansion ap- 
proach [4]. The Coulomb-Sturmian (CS) functions are defined by 



(r\nl) = 



1 1/2 



(n + 2/ + l)! 



(26r) z+1 exp(-6r)L^ +1 (26r), 



21+1, 



(18) 



with n and / being the radial and orbital angular momentum quantum num- 
bers, respectively, and b is the size parameter of the basis. The CS functions 
{|n/)} form a biorthonormal discrete basis in the radial two-body Hilbert 
space; the biorthogonal partner is defined by (r\nl) = (r\nl)/r. Since the 
three-body Hilbert space is a direct product of two-body Hilbert spaces an 
appropriate basis can be defined as the angular momentum coupled direct 
product of the two-body bases 

\nul\)i = \nl)i <g) \v\)i, (n,u = 0, 1,2, . . .), (19) 

where \nl)i and are associated with the coordinates Xi and yi, respec- 

tively. With this basis the completeness relation takes the form (with angular 
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momentum summation implicitly included) 

N 

1 lim V |nWA)i i(ra4A| = lim if. (20) 

n,u=0 



Note that similar bases can be constructed for fragmentations 2 and 3 as well. 
We make the following approximation on Eq. (16) 

\1>) = G?l?v[ a) lfPl?\1>), (21) 



(s) 

i.e. the operator v\ pV in the three-body Hilbert space is approximated by a 
separable form, viz. 



N 



v\ s 'pV= lim l^v[ s ' P Vl^ 

N—hx> 



]T |nWA)i ^ i(nVZ'A'|, (22) 

ri,i/,n',i/=0 

where uf^ = i(nul\\v[ s ^pV\n'u'l'X')i. The compactness of the equation and the 
completeness of the basis guarantee the convergence of the method. Utilizing 
the properties of the exchange operator V these matrix elements can be written 
in the form = p x (— ) 1 ' i{nvl\\v^\n'v'l'\')2- 

With this approximation, the solution of Eq. (16) turns into solution of matrix 
equations for the component vector ip = i(ni4\\i/;) 

{^ ,) W]- 1 -eS% = > (23) 

where G? = i(^|Gf Vz//'A')i. A unique solution exists if and only if 

D(z) = det{|G? ) (z)]- 1 - v[ s) } = 0. (24) 

The Green's operator Grp is related to the Hamiltonian H[ l \ which is still 
a three-body Coulomb Hamiltonian, seems to be as complicated as H itself. 
However, has only 1-type asymptotic channels, with asymptotic Hamil- 
tonian 

H 1 = H° + vf. (25) 
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Therefore, in the spirit of the three-potential formalism [5], G[ can be linked 
to the matrix elements of G\(z) = (z — Hi)' 1 via solution of a Lippmann- 
Schwinger equation, 

(Si )" 1 = (Si) -1 "Ei, (26) 

where G± = l (^^\G l \rM\') l and U_y = ^nvlX^ + )|nVZ'A')i. 

The most crucial point in this procedure is the calculation of the matrix el- 
ements G_i, since the potential matrix elements and \J_y can always be 
evaluated numerically by making use of the transformation of Jacobi coordi- 
nates [6]. The Green's operator G\ is a resolvent of the sum of two commuting 
Hamiltonians, H\ = h Xl + h yi , where h Xl = h Xl + ff(xi) and h yi = h yi , 
which act in different two-body Hilbert spaces. Thus, G\ can be given by a 
convolution integral of two-body Green's operators, i.e. 

Gi(z) = ^-fdz'g Xl (z - z>) g yi (z'), (27) 
c 

where g xi (z) = ( z ~ ^xj -1 an d 9yi( z ) — ( z ~ The contour C should 

be taken counterclockwise around the continuous spectrum of h yi such a way 
that g Xl is analytic on the domain encircled by C. 

In time-independent scattering theory the Green's operator has a branch-cut 
singularity at scattering energies. In our formalism G\(E) should be under- 
stood as G\{E) = \im e ^oGi(E + ie), with e > 0, and E < 0, since in this 
work we are considering resonances below the three-body breakup threshold. 
To calculate resonant states G\(E + ie) has to be continued analytically to 
negative e values. Before doing that let us examine G\(E + ie) with e > 0. 
Now, the spectra of g yi {z') and g xl {E + ie — z') are well separated and the 
spectrum of g yi can easily be encircled such that the spectrum of g Xl does not 
penetrate into the encircled area (Fig. 3). Next, the contour C is deformed 
analytically in such a way that the upper part descends into the unphysical 
Riemann sheet of g yi , while the lower part of C is detoured away from the cut 
(Fig. 4). The contour still encircles the branch cut singularity of g yi , but in 
the e — > limit it now avoids the singularities of g xi . Moreover, by continuing 
to negative values of e, in order that we can calculate resonances, the pole 
singularities of g xi move onto the second Riemann sheet of g yi (Fig. 5). Thus, 
the mathematical conditions for the contour integral representation of G±(z) 
in Eq. (27) can be fulfilled also for complex energies with negative imaginary 
parts. In this respect there is only a gradual difference between the bound- 
and resonant-state calculations. Now, the matrix elements G a can be cast in 
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g. (Kii?-?') 

!i \e=. 


■ 

z' 




g y Jz') 



Fig. 3. Analytic structure of g xl (E + ie — z') g yi (z') as a function of z' , e > 0. The 
Green's operator g yi (2') has a branch-cut on the [0, oo) interval, while g Xl (E+ie—z 1 ) 
has a branch-cut on the (— oo, E-\-ie] interval and infinitely many poles accumulated 
at E + ie (denoted by dots) . The contour C encircles the branch-cut of g yi . In the 
e — > limit the singularities of g Xl (E + ie — z 1 ) would penetrate into the area covered 
byC. 




; 



Fig. 4. The contour of Fig. 3 is deformed analytically such that a part of it goes on the 
unphysical Riemann-sheet of g yi (drawn by broken line) and the other part detoured 
away from the cut. Now, the contour avoids the singularities of g Xl (E + ie — z') even 
in the e — > limit. 

the form 

^ = hf dz ' 9 -^ z - z ' )g -y^ (28) 

c 



where the corresponding CS matrix elements of the two-body Green's opera- 
tors in the integrand are known analytically for all complex energies (see [5] 
and references therein). It is also evident that all the thresholds, correspond- 
ing to the poles of g , are at the right location and therefore this method is 
especially suited to study near-threshold resonances. 
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Fig. 5. In the e < case, what is needed to calculate resonances, the poles of g Xl 
which lie above the branch point of g yi go to the unphysical sheet of g yi (denoted 
by shaded points), while the others remain on the physical one. 

5 S'-wave resonances in positronium ion 



To calculate resonances we have to find the complex zeros of the Fredholm 
determinant (24). To locate them we calculate D(z) along the real energy axis 
with small step size. In the vicinity of zeros the D(z) exhibits violent changes. 
Then we have a good starting point for the zero search. This way we can easily 
find resonances, at least the narrow ones. To find broad resonances one should 
proceed by calculating D(E — ie) with finite e > 0. 

For the parameters of the splitting function we take y° = 50ao and v = 2.1, a® 
is the Bohr radius. To select out the physical solutions we vary the short- and 
long-range potentials by taking x° = 18a , 25a and 30a , respectively. We can 
also vary the potentials by adding and subtracting a short-range term 



v{ s) = v^C(xi,yi)+v C(x l7yi ) (29) 
v? = v? [1 - C(zi, yx)} - v C(xi, yi), (30) 

respectively, while keeping a fixed x° = 18ao value. This new kind of split- 
ting goes beyond the Merkuriev's original suggestion, but since the charac- 
ter of v( s > and remained the same all the nice properties of the origi- 
nal Faddeev-Merkuriev equations are retained. In these calculations we used 
v = 0,0.01,0.02 values. 

We found that some solutions, especially their widths, are very sensitive to the 
change of either x° or vq parameters. We regard them spurious solutions. Those 
solutions, given in Table I, were stable against all changes of parameters and 
we can consider them as physical resonances. We recovered all the resonances 
presented in Ref. [3] with very good agreements, but besides that we found 10 
more resonances, all of them are in the vicinity of some thresholds. 
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State 


Ref. [3] 


Our work 


1 ce 


T? 

— Fj r 


1 


171 

— Cj t 


1 




0.1520b08 


0.00008b 


0.15192 


0.000042b 




U.IzyJU 


U.UUUU2 


0.12727 


0.0000084 








n 1 on 

0.1251 


C\ AAAAAO 

0.000002 




U.U7l)00<3 


U.UUU15 


0.70DDD 


0.000074 




U.Uoyoy 


(J.UUU11 


0.059b82 


0.000052b 








0.0504 


0.00003 




0.04045 


0.00024 


0.4042b 


r\ nnn 1 on /i 

0.0001294 




n no £n 
0.0350 


n nnno 
0.0003 


0.035020b 


n nnm o 
0.00013 




0.03403 


0.00034 


0.0346234 


0.000158b 








0.032b3 


n nnn i 

0.0001 








n n^i 


n nnnn7 

U.UUUU ( 




U.UZDo 


n nnn/i ^ 

U.UUU4D 


U.UZOUD 


n nnm n/i 

U.UUU1U4 




U.UZo4o 


n nnm a 

U.UUU14 




n nnnn/i 

U.UUUU4O0 


3oe 

kJ 


771 

— _C/ r 


■p 
l 


771 

— Hi r 


l 




0.1270b 


0.00001 


0.127050 


0.0000000028 








n 1 on 

U.lzol 


0.000000001 




0.05873 


0.00002 


0.05872 


n, AAAnnAi oro 

0.0000001852 








0.05bl 


0.0000002 








n n^^ 


n nnm 

U.UUUl 




0.03415 


0.00002 


0.0342018 


0.00000070 








0.03237 


0.0000011 








0.03172 


0.0000002 








0.031035 


0.0000938 



Table 1 

L = resonances of Ps~. The energies and widths are expressed in Rydbergs. 
6 Conclusions 

In this article we pointed out that in the Faddeev-Merkuriev integral equation 
approach of the three-body Coulomb problem the complex-energy spectrum 
is contaminated with spurious solutions. The spurious solutions, however, are 
sensitive to the splitting of the potential into short- and long-range terms. 
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This offers an easy way to select the physical solutions. 

We solved the integral equations by using the Coulomb-Sturmian separable 
expansion approach. This method gives an exact description of the thresh- 
old phenomena, thus the method is ideal for studying close-to-threshold reso- 
nances. In the e _ e + e~ system we located 10 new resonances, all of them are 
in the close vicinity of some threshold. 

This work has been supported by the NSF Grant No.Phy-0088936 and OTKA 
Grants under Contracts No. T026233 and No. T029003. 
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